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J/3 , Abstract 

In this paper, a multiplicity preserving triangular set decomposition algorithm is 
proposed for a system of two polynomials. The algorithm decomposes the variety de- 
fined by the polynomial system into unmixed components represented by triangular sets, 
which may have negative multiplicities. In the bivariate case, we give a complete algo- 
rithm to decompose the system into multiplicity preserving triangular sets with positive 
multiplicities. We also analyze the complexity of the algorithm in the bivariate case. 
We implement our algorithm and show the effectiveness of the method with extensive 
experiments. 

Keywords. Triangular set decomposition, multiplicity preserving decomposition, ex- 
tended Euclidean algorithm. 



X. 

1 Introduction 



Decomposing a polynomial system into triangular sets is a classical method to solve polyno- 
mial systems. The method was first introduced by Ritt [2UJ and revised by Wu in his work 
of elementary geometry theorem proving [24J [25]. There exist many work about this topic 
[DI3llligiSl[naiIDIiaiI3II3iniIISll22l|2Zl. The main tool to decompose a polynomial 
system is pseudo-division. In most existing triangular decomposition methods based on the 
pseudo-division algorithm, one need to deal with the initial of certain polynomial(s), say h, 
which will bring extraneous zeros. Usually, one decomposes here the system into two systems 
corresponding to the cases h = and h ^ 0. Doing so, the number of the systems increases 
quickly. Moreover, this leads to some repeated computations which can be avoided. 

Another reason why we consider the topic is that the multiplicity of a component or a 
zero of a polynomial system is an important information which helps us to obtain a further 
understanding of the structure of the variety defined by the polynomial system. 
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Most triangular set decomposition algorithms do not preserve the multiplicities of the 
zeros or the components. One approach to remedy this drawback is to decompose the 
polynomial system into triangular sets first and then recover the multiplicities. Li proposed 
a method to compute the multiplicities of zeros of a zero-dimensional polynomial system 
after obtaining a triangular decomposition of the system [12]. Recently, Li, Xia, and Zhang 
proved that the characteristic sets in Wu's sense for zero-dimensional polynomial system is 
actually multiplicity preserving with a minor modification [TJ] ■ They also gave a multiplicity 
preserving decomposition, but some of the components are not in triangular form. 

In this paper, we use the concept of multiplicative variety, that is, the components and 
their multiplicities in the original polynomial system. We consider not only the components 
themselves but the multiplicities of these components. During the decomposition, the initials 
bring some extraneous multiplicative varieties in each pseudo-division step. We record them 
during the computation and remove them later, which helps us to recover the multiplicative 
varieties of the original system. We also avoid some repeated computation during the decom- 
position. Currently, the theory is complete for polynomial system with two polynomials. In 
particular, we provide a method to compute the multiplicative-zeros of a zero-dimensional 
bivariate system with two polynomials. We also analyze the complexity of the algorithm 
under some conditions. 

Kalkbrener's method for zero-dimensional bivariate polynomial system is similar to our 
method [17] . But his method is not multiplicity preserving. And our method is in a different 
sense: we remove the extraneous zeros from the system. 

The paper is organized as below. In the next section, we provide some properties of prim- 
itive polynomial remainder sequences. In Section 3, we provide the theories to decompose a 
polynomial system with two polynomials into triangular sets which preserve the multiplicities 
of the components of the original system. We provide a multiplicity preserving algorithm to 
decompose a zero-dimensional bivariate polynomial system into triangular sets in Section 4. 
The complexity of the algorithm under some conditions are analyzed. Algorithms and ex- 
amples are used to illustrate the effectiveness and efficiency of our method. We also compare 
our method with other related methods. We draw a conclusion in the last section. 

2 Primitive Polynomial Remainder Sequence 

In this section, we introduce some basic properties for primitive polynomial remainder se- 
quences. In fact, there are many references for this topic, in particular [2|I161[TT]. We modify 
the procedure for our own purpose. 

Let K be a computable field with characteristic zero, such as the field of rational numbers 
and K[yi, . . . , y n ] the polynomial ring in the indeterminates y%, . . . , y n . 
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Let p G K[xi, . . . , x n , x]. We define 

Cont(p, x) = gcd(coeff(p, x l ), i = 0, 1, . . . , deg(p, x)), 
Prim(p, x) = p/Cont(p, x), 

where coeff(p, x' 1 ) means the coefficient of x % in p and deg(p, x) means the degree of p in x. 
p is called primitive w.r.t. x if Cont(p, x) = 1. 

The pseudo-division can be extended to the following form. 

Lemma 2.1 Let f,g G K[x\, . . . ,x n ,x], deg(/,x) = d±, deg(g,x) = d 2 , d\ > d 2 , and 
gcd(/, g) = 1. There exist q,r £ K[x±, . . . , x n , x] such that 

l 5+1 f + qg = r, (1) 

where I is the leading coefficient of g in x, 5 = d\ — c?2, deg(g, x) > deg(r, x). Furthermore, 
q has the form: 

q = ltx + s, (2) 

where t G K[x±, . . . , x n , x], s G K[x±, . . . , x n ]. Moreover, if r\ = Cont(/, x),r2 = Cont(<7, x), 
then 

n\q, r^ d2 \q, ri \r,r^ d2+1 \r. (3) 
Proof. Write f,g as univariate polynomials in x, 

/ = ai x dl +a 2 x <il - 1 + --- + a dl+ i, 
g = b lX d2 +b 2 x d2 ~ 1 + --- + b d2+1 . 

To eliminate the terms of / with degree d\ in x, we have 

Tq(x) = b\ f + qo g = ho x dl ~ x + lower powers in x, 

where qo = —aix dl ~ d2 ,ho = b\ a 2 — a\b 2 . It is clear that r\\qo, since r\\a\, r®(= l)|(?o and 
ri\To,r2\To. So the lemma holds when 5 = 0. Note that r = To when 5 = 0. Now, we need 
to eliminate ho * x dl ~ x from Tq{x). If ho / 0, 

Tx(x) = b 1 T (x)-(b 1 a 2 -a 1 b 2 )x d ^- 1 g 

= b\f + fa g - (6i a 2 - ai b 2 ) x^'^g 

= b\f + q ig 

= h\ x dl ~ 2 + lower powers in x, 

where h\ G K[x±, . . . , x n ]. Each term of q\ contains a factor of the form a« bj. So n|gi, r2\qi 
and ri|Ti,r||Ti. And <?i = -&iaix dl ~ d2 - (&i a 2 - aib 2 ) x dl_d2_1 . If /i = 0, the results is 
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still true. So the lemma holds when 5 = 1. Assuming that the lemma holds for the cases 
5 < i, then we have 

Tj(x) = b{ +1 f + qj g = hj x dl ^^~ l + lower powers in x, 

Note that deg((?j_i, x) > deg(gj,a;) and the lowest power of q<j—i in x is larger than d\ — d 2 —j 
and r%\hj, r^ +1 \hj. Then 

Ti + \{x) = b l ^~ 2 f + qi + i g = hi + \ x dl ~ l ~ 2 + lower powers in x, 

We can similarly derive qi+i = b\qi — hi x dl ^ d2 ^ 1 ^ 1 . When 5 = i+1, we have Oj+i = b\qi — hi. 
So <7j+i has form ([2]) since the lowest power of qi is larger than 0. And n |oj_|_i since T\\qi,r\\hi. 
So r\\Ti+\{x). r l 2 +1 \qi+i since r 2 \b\,r 2 \qi and r 2 +1 \hi. So r l 2 +2 \Ti + i(x). So the lemma holds 
for 5 = i + The lemma is proved. | 

Corollary 2.2 Let f,g € K[xi, . . . ,x n ,x] be primitive, di = deg(/, x) > d 2 = deg(g,x), 
and gcd(/, g) = 1. Regard f,g as univariate polynomials in x. Then there exist an m £ 
. . . ,x n ] suc/j i/iai 

™/ = <?5 + r, 
and gcd(m, g) = gcd(m,r) = 1. Furthermore, 

(mf,g) = (g,r), 

where (P) represents the ideal generated by P. 

Proof. The corollary is obvious. | 

Corollary 2.3 Cont(<7, x) = 1 if gcd(l,s) = 1 and d\ > d 2 , where I and s are from fTp and 

Proof. Regard /, g as univariate polynomials in x, and q a polynomial in x and dj, bj, where 
i = 1, . . . , di + 1, j = 1, . . . , d 2 + 1. Let r 2 \g and r 2 £ K[x±, . . . , x n \. From Lemma I2TTI r 2 \q 
if di > c?2 • So r 2 \s. Since rs r2|gcd(Z,s). We have r 2 = 1 if gcd(/,s) = 1. The lemma is 
proved. | 

The above result is a necessary condition to check whether g has factors in J^C[xi, . . . , 

Lemma 2.4 Let fi,f 2 S K[xi, . . . ,x n ,x], d\ = deg(/i,x) > d 2 = deg(f 2 ,x). Assume that 
Cont(/j,x) = l,i = 1,2. Applying the extended Euclidean algorithm for fx, f 2 w.r.t. the 
variable x, we obtain a polynomial sequence {fi,f 2 , ■ ■ ., fk+2} such that 

mfi + qifi+i = m i ^ 1 pif i+2 ,i = 1,. . . ,k, (4) 



4 



where m = l,p k = l,m i} pi, fk+2 € K[xi, . . . , x n ], q t G K [xi, . . . , x n , x], i = l,...,k, and 
Cont(/j,x) = 1(1 < i < k + 1), gcd(mi,pi) = 1. 

Proof. We prove the lemma by induction on i. When i = 1, from Lemma l2.1( there exist 
q G K[x\, . . . , x n ,x], r G K[x±, . . . , x n ,x] such that l 2 +1 fi + qH = r, where U is the leading 
coefficient of fi in x, 5 = d\ — d 2 - Let t = gcd(/2 +1 , q), mi = , q\ = |. Let p\ = Cont f ( r '' :E ) 
and /3 = Prim(r, x). It is clear that gcd(mi,pi) = 1- Assume that for 1 < i < i, dH) holds. 
Denote di = deg(/j). For j = i, we have lf+lfi + ?t/i+i = r i+2 , where 9 = di - d i+1 . If 
wjj-i is a factor of rj+ 2) set p\ as the product of all the factors of in if [xi, . . . , x n 1. Let 

/i = gcd(Z i + 1 1 ,p.)- Then m i = ft = %>Pi = gcd(mi,pi) = 1. If mj_i is not a factor 
of rj + 2, we can multiply g = gc d(m *~* r +2 ) *° ^ e * wo s ^ es °f the equation. Then doing the 
same operation as before, we can derive rriifi + qifi+i = rrii-iPifi+2 which satisfies all the 
conditions. This proves the lemma. I 

Remark: In most cases, we have gcd(mj,(/j) = 1 and pj = 1 which helps us to design 
efficient algorithms. 

Corollary 2.5 Let f u f 2 G K[xi, . . . , x n , x], gcd(/i, / 2 ) = 1, and Cont(/j,x) = l,i = 1,2. 
From the extended Euclidean algorithm, we can obtain 

mifi + ft/i+i = i = 1, • • • , (5) 

(mfiJi+i) = (fi+i,9i fi+2), (6) 

w/iere mi,£i G K[xi, . . . , x n ], gcd(mi,gi) = I, g k = 1, gcd(m fc , g k f k+2 ) = 1, and / i+2 (l < 
i < — 1) are primitive. 

Proof. From Lemma 12.41 we have (HJ). Note that gcd(mj,j>j) = 1. Let h = gcd(mj, mj_i), 
denote m; = = = mi ~ m . Then we have gcd(mj,gj) = 1. Since / fc+2 G 

if [xj, . . . , x n ], we can set g k = 1. We can delete gcd(m k , g k f k+ 2) if it exists. ([6]) is ob- 
vious. So the corollary holds. I 
The following corollary is clear and useful. 

Corollary 2.6 We can rewrite {5|j and f^j as below. 

rrhfi + qifi+i = —^Pifi+2,i = l,...,k, (7) 

Wi 

(mifi, fi+i) = [f i+ i,— ! —pifi + 2), (8) 
where Wi is a factor ofrrii-i, gi = 1711-1 Pi, and p k = 1. 
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3 Triangular Decomposition of Two Polynomials 



In this section, we will give the method to decompose a system of two polynomials into 
triangular sets. We need the concept of multiplicity variety. 

Definition 3.1 ([9] pp. 129-130) Let be an unmixed variety of dimension d in a projective 
space S n of dimension n over K . And 

1=1 

where V\ is an irreducible variety of dimension d and order g^. Let Fi(iiQ, . . . ,Ud) be the 
Chow form (see J2}/ pp.32) of vj ; which is irreducible over K and of degree gi in uj = 
(ujo, u jn ), for j = 0, . . . , d. The form 

F{u ,...,u d ) = II , * =1 [F l (uo,...,u d )} a \ (9) 

where ax, . . . , are positive integers, satisfies the conditions for a Chow form of an algebraic 
variety which, regarded as a set of points, coincides with V4. We consider a new entity, 
consisting of the variety Vd associated with the form F{uq, . . . , Ud), for a given choice of the 
exponents ax, ■ ■ ■ ,Oh, denoted as M^. We write 

h 

M d = ^o i M«, (10) 
i=i 

(i) U) 

where M\ corresponds to V) . We call a multiplicative variety and the multi- 

(i) 

plicity o/M\ . Especially, we call a multiplicative-zero set when d = 0. 

Remark: Since an affine variety can be easily transformed into a projective variety, we will 
consider directly affine multiplicative varieties in K n in this paper. And we assume that the 
system has no solutions at 00. 

Theorem 3.2 (J3j[ ppl60) M a , Ml, Mi are unmixed multiplicative varieties with dimension 
a,b respectively. The intersection o/M a and M^ + M^' are with dimension a + b — n. So are 
the intersections o/M a and Mi, M a and M^'. Then we have 

M a ■ (M' b + M'l) = M a ■ M' b + M a • M?, 

where ■ represents the intersection of two multiplicative varieties, which preserves the multi- 
plicities of each intersection component (for more details see \9j, ppl58-160). 



6 



Definition 3.3 (f^j pp.139) Let I be a zero dimensional ideal in K[x±, . . . ,x n ] such that 
the variety V(I) defined by I consists of finitely many points in K , where K is algebraic 
closure of K , and assume p = (01, . . . , a n ) G V(I). Then the multiplicity of p as a zero of 
I, denoted by m(p), is the dimension of the ring obtained by localizing K[x\, . . . ,x n ] at the 
maximal ideal M = I(p) = {x\ — a\, . . . , x n — a n ) corresponding to p, that is: 

m(p) = dim-^^i, . . . ,x n ] M /IK [x 1} . . .,x n ] M . 

Lemma 3.4 (J^ pp. 144 ) Let {ft,..., /„} C K[ x\, . . . , x n ] be zero- dimensional, and have 
total degrees at most d\, . . . , d n and no solutions at 00. If fo = uq + u\ x\ + . . . + u n x n , where 
uq, . . . ,u n are independent variables, then there is a nonzero constant C such that the Chow 
form of {fx,.. . ,/ n ) is 

Res Mlv .. irfn (/o,...,/ n ) 

where p = (£1, . . . ,£ n ) G • • • , f n ) and m(p) is the multiplicity of p in (f x , ... , f n ). 

The lemma illustrates the relationship between Chow form and the multiplicity of a 
point of a zero-dimensional polynomial system. The lemma tells us that m(p) is the multi- 
plicity of the corresponding irreducible zero-dimensional component of the zero-dimensional 
polynomial system {f\, . . . , f n } when d = in Definition 13.11 

From Theorem 13.21 we have the corollary below. 

Corollary 3.5 Using the notations as Corollary \2.S\ we have 

M(mf,g)=M(m,g)+M(f,g). 

The following lemma is important to our algorithm. 

Lemma 3.6 Let f,g,r,m be as Corollaru \2.$\ We have 

M(f,g)=M(g,r)-M(m,g). (11) 

Proof. By Corollary 12.21 (mf,g) = (g,r), so we have 

M(mf,g)=M(g,r). (12) 

We can find that m f,g define a multiplicative variety with dimension n — 2 since 
gcd(m /, g) = 1. Note that /, g both are primitive and m G K[x%, . . . , x n \. And /, g and m, g 
both are n — 2 dimensional. 
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We are going to prove that 



M(f,g)=M(mf,g)-M(m,g). (13) 

Assume that the Chow form To of M.(mf,g) is as Q and the order of Wl(mf,g) is 
Go = Yli=i a i 9i- From (|3.5|) . we can assume that the Chow form of M(m, g), M(/, g) are 
T\ = n-i =1 [Fi(u , . . . , u d )] a > and T 2 = U^ k+1 [Fi(u , u d )] a \ respectively, where 1 < k < h 
and d = n — 2. And the orders are G\ = J2i=i a i 9i an< ^ G 2 = J2i=k+i a * ft> respectively. We 
define the Chow form T' 2 of M(m /, <?) — M(m, <?) as Tq/T\ and the order as G2 = Go — Gi. 
Since the Chow form of the algebraic variety is not equal to zero, the definition is well 
defined. Here (m,g) defines some components of (mf,g), including the multiplicities of 
the components. Thus M.(mf,g) — M,(m,g) still has positive exponent for each simplified 
component. And we have 

T> 2 = T$ =l [F i {u Q ,...,u d )P/{Tl k i=l [F^ 

h k h 

G' 2 = S ^aigi-^ j a i g i = aigi = G 2 . 

i=l i=l i=k+l 

So (fT3j) holds. Combining (fT2|) and (fTB"]) . we have (jTTj) . This ends the proof. | 

Definition 3.7 A multiplicity preserving triangular decomposition of a polynomial 
system £ is a group of triangular sets {T^,T^ ,i = 1, . . . ,m + ,j = 1, . . . ,m~} in multiplica- 
tive variety sense such that 

M(E) = J^M(T+) -J^M(7T). (14) 
i=i j=i 

Remark: We will show that a multiplicity preserving triangular sets decomposition exists 
for systems with two polynomials in the rest of the paper. Note that for a zero-dimensional 
polynomial system, the existence of f|14[) is obvious. The existence of (I14p for for general 
case (dimension mixed, more polynomials) is our future work. 
The following is a key result of the paper. 

Theorem 3.8 Let fi, f 2 € K[x±, . . . ,x n ,x] such that gcd (fi, f 2 ) = 1 and Cont(/j, x) = 
l,i = 1,2. Then 

k k 

M(/i, h) = M(/ fc+ i, A+a) + £ M( 5i , - £ M(mi, (15) 

8=1 1=1 

where fi,gi,rrii are defined in Corollary \2.5[ 
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Proof. From ([6]), we have 

M{m i f i ,f i+l ) = M{f i+l ,g i f i+2 ). 

Then by Corollary 13,51 f° r 1 < i < A;, we have 

M(m i5 f i+l ) + M(/i, = Mfo, + M(/i+i, / i+2 ), (16) 
M(/i, fi+i) = M(/ m , / i+2 ) + M( 9i , f i+1 ) - M(m«, (17) 

So we have 

M(/!,/ 2 ) = M(f 2 ,f 3 )+M(g 1 ,f 2 )-M(m 1 J 2 ) 

= M(/ 3 , / 4 ) + MG/i, / 2 ) + M( 52 , f 3 ) - M(mi, / 2 ) - M(m 2 , / 3 ) 

= M(/ fe+1 , / fc+3 ) + £ M(flj, - M K. 

i=l i=l 

Remark: The decomposition is about the (n — 2)-dimensional component of M(/i,/ 2 ). 

Corollary 3.9 Use the notations as Corollary 1 2. 6\ we have 

k 



M(/i, f 2 ) = M(/ fc+1 , A +2 ) + X)(M(p(, /i+i) - M(wi, f i+l )) 

i=l 

fe-i 

- YVM(m i; ^i) + M(mi,pi) - M(rm, ?*)) - M(ro fc , / fc +i). (18) 

Proof. From Corollary I2.6| we have u>i = mi-ipi. So by (|12p and Corollary 13.51 we have 

M(<7i, /j+i) = M(pi, + M(mi_i, - M(m ( , (19) 
This helps us simplifying the computation. By (|17p . we have 

, / a ) = MCfk+i, fk+2) + £(M(pi, + M(m4_i, /i+x) - M(m ( , - ]T M(m 4 , 

i=l i=l 
fe fc fe-1 

= M(/fc+i, /jfc+2) + M (Pi>fi+i) - E M ( w " + M ( m °' /2 ) + E M(m i5 / i+2 ) 

i=l i=l i=l 

fe-1 

- ]T M(m<, /i+i) - M(m k ,f k+1 ). (20) 

i=l 

, . , rrii-i , 

mfi + qiji+i = Piji+2, 

Wi 



From 



we have 

777 ' 1 

M( mi , qi f l+1 ) = M{m h -^± Pi f l+2 ). 

Wi 

By (fT2|) and Corollary 13.51 we have 

777 ' 1 

M(mi,/ i+ i) = M(m i ,f i+2 )+M(rn i ,^^)+M(rn i ,p i )-M(rn i ,q i ). 

Wi 

And m = 1, so M(m , A) = 0- Then we have (flB]). | 
Remark: The components M(mj, ^T 1 ), M(mj,pj)), and M(mj, only involve polynomials 
in K[x\, . . . ,x n }. Note that by Lemma 12.41 the coefficient of % in x* for i > is zero when 
rrii = 0. These components can also be decomposed into triangular sets recursively. This 
corollary is very important since it provides a method to eliminate the main variable x in 
fas, which simplifies the decomposition. We can obtain another interesting phenomenon 
from simple observation, that is, the degree of all the resulting polynomials is bounded by 
the square of the degree of fa, fa. 

Example 3.10 Consider the system [fa, fa] = [x 2 + y 2 + z 3 — 1, xz 2 — zy + 1] with Corollary 
\3.9\ under the variable order x -< y ~< z, 

t 4,2 2 2 , 2 

J3 = x+xy— x— xz — y + yz, 

fa = 1 - 3x 3 y - 3xy 3 + 3xy + x 2 y 3 + y 5 - y 3 + x 7 + 2x 5 y 2 - 2x 5 + x 3 y 4 - 2x 3 y 2 + x 3 . 

mi = x 2 ,qi = xz + y, tuq = wi = p\ = 1, 7712 = (— x + y 2 ) 2 ,q2 = —x 2 z + xy 2 z + 2 xy — y 3 — 
x 5 — x 3 y 2 + x 3 . W2 = P2 = 1- By Corollary \3.fA we have the following decomposition. 

M(/i, fa) = M(/ 4 , fa) + M(mi, q x ) - M(m 2 , / 3 ), 

where M(m 1; gi) = 2M(x,y) and M(m 2 ,/ 3 ) = 2M{x,y) + 2M{x - l,y - 1) + 2M(/ii,/i 2 ), 
where h\ = x e + 3x 5 + 2x 4 + x 2 + x + 1, /12 = y — x 4 — x 3 + x 2 . T/ius we /iaue 

M(/i, / 2 ) = M(/ 4 , / 3 ) - 2M(x - 1, y - 1) - 2M(/»i, 

iVoie i/iai the component with negative multiplicity cannot be removed if using triangular 
form. 

4 Multiplicity Preserving Decomposition for System of Two 
Bivariate Polynomials 

In this section, we will consider the triangular decomposition of a zero-dimensional bivari- 
ate polynomial system with two polynomials, that is, S = {/, g} C K[x,y]. The method 
provided here is complete for a zero dimensional bivariate polynomial system with two poly- 
nomials. When E is zero dimensional, M(S) defines a multiplicative-zero set. 
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4.1 Algorithm 

Lemma 4.1 Using the similar notations as Corollary \2.(k «/{/i,/ 2 } is zero- dimensional, 
we have 

k k 

M(/i,/ 2 ) = M(/ w ,/w) + £M( Pi) / i+1 ) - ^M^,/^) - MK,/ W ). (21) 

i=l i=l 

Proof. The lemma is a consequence of Corollary 13.91 Note that M(m.j, ) = M(mj,j?j) = 
M(mj,gj) = since gcd(mj, pi), gcd(m,j, rr ^~ 1 ), gcd(mj,qj) are constant. 
The following corollary is useful. 

Corollary 4.2 If Wi{\ < i < k) are constants and fk+i = h(x)y t + lo(x) for t > and 
Zo,'i € K[x] in nn\). we have 

k 

M(/i,/ 2 )=M(/ Jfc+lj / fc+2 ) + ^M(pj,/ i+1 ). (22) 

i=l 

Furthermore, if Pi(l < i < fc) are constant, we have 

M(/i,/ 2 )=M(/ fc+ i,/ fc+2 ). (23) 

Proof. 

Since = h(x)y t + / (a;) and Cont(/ fc+ i, y) = 1, M(m fc ,/ fc+ i) = 0. Note that m fe is 
a factor of Z" for some positive integer n. Thus M(mo,/2) = 5 M(mjt, fk+i) = 0- So from 
([H]), we have ((25]). If = 1, ((23]) is a consequence of (|55J. | 

Now we will consider the complexity of our method under the condition of Corollary 14.21 
We consider this case because it is usually the case for almost all zero-dimensional bivariate 
polynomial system with two polynomials. So the result is interesting. At first, we need 
to introduce some notations, which can be found in [7]. Let C(f) bound the bitsize of the 
coefficients of / € K[x,y] (including a bit for the sign). We assume lg(deg(/)) = 0{C(f)). 
For a € Q, C(a) is the maximum bitsize of a's numerator and denominator. Let M(r) 
denote the bit complexity of multiplying two integers of size r, and M(d, r) the complexity 
of multiplying two univariate polynomials of degrees < d and coefficient bitsize < r. Using 
FFT, M(r) = O b (t) and M(d,r) = d B (dr). 

Lemma 4.3 $[W$ Let f,g G Z[x], deg(/),deg(y) < n, and C(f),C(g) < r. We can 
compute gcd(/, g) in Ob(ji 2 t). 

Lemma 4.4 P7[ {T^ Let f,g G Z[x,y], deg(/), deg(p) < n, and C(f),C(g) < r. We can 
compute the subresultant sequence of f and g in Osin^r). 
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Theorem 4.5 If Wi(l < i < k) are constant and fk+i = h(x)y t + Iq[x) for some positive 
integer t, Iq,1\ £ K[ x ] ^ n (E2P; an d K = r L, we can decompose a zero- dimensional bivariate 
system with two polynomials into multiplicity preserving triangular sets in OB(n 7 r). 

Proof. We can compute a subresultant sequence of / and g at first. It can be computed 
in 0e(n 6 r) by Lemma 14.41 Then we simplify each pseudo-division step to derive (J7J) from 
the highest degree of the sequence in y to the lowest degree. Let {iq, . . . , ^+2} be the 
subresultant sequence of / and g. We need only consider the case of regular subresultant 
sequence since the complexity of the regular case also bounds the degenerate case. Consider 
the formula 



Assume that we have computed the contents of Fi and -Fj+i, say rj,rj + i. F\,F<z are /i,/2- 
And the contents of / or g can be computed in Os(n 3 r) by Lemma [4. 3 \ which can be ignored 
comparing to 0B(n 6 7~). For each Fi, it is well known that deg(-Fj) < n 2 ,C(Fi) = 0(nr) (for 
reference see [7]). And deg(Fj, y) < n— 1 for i > 3. Thus for any coefficient of Fi, say h G Z[x], 
we have deg(h) < n 2 ,C(h) = 0{nr). So to compute the content of Fi with deg(Fj,y) = h, 
we need to compute at most h gcd each in Os(n 5 r). Let r-i + 2 = Cont(Fj + 2, y). In order 
to derive ((7J), we need to delete gcd(l? +1 Fi, ifF^) from the two side of ([24"]) . So we 

need to bound gcd(/? +1 Vi, if ^+2). Note that deg(s) < n 2 ,C(s) = 0{nr) holds for s = h or 
s = li + \ and we can not optimize the degree of h+i- But we can compute r = gcd(/j, k+i), 
which is bounded by Ob(h 5 t). And gcd((^-) 2 , ri + 2) can be bounded by 2 0s(n 5 r) as below. 
We can compute w = gcd((^-),rj+2), and then gcd((^), gcd((^-) 2 , rj) is also bounded 

by 2 d B (n 5 r). So to obtain © from (j2U) for deg(F i5 y) = fc, we need (fc + 5)O e (n 5 r). Then 
we can decide rrii,gi in ([5]) by two divisions. Since tUj's are constant, Thus, we can 

obtain pj by one division. When /i changes from n to 1, we can bound it by n " 2 " 9n C ) g(?7- 5 r), 
that is, Ob{it. 7 t). Then the total complexity is Os(n 7 r). | 
Remark: For many f,g G the last two elements of the subresultant sequence 

-Ffc+i, form a multiplicity preserving triangular decomposition of f,g. Thus, we can 
compute the decomposition in Ob(b 6 t). 

The lemma gives a multiplicity preserving triangular decomposition of a bivariate polyno- 
mial system. But there exist some triangular sets with negative multiplicities. The following 
results gives an algorithm to remove the triangular sets with negative multiplicities. 

Theorem 4.6 There exists an algorithm to decompose a zero- dimensional bivariate polyno- 
mial system {/i,/2} C K[x,y] into a set of triangular sets, such that 

N 



if+iFi + QiFi+i — l 2 Fi + 2- 



(24) 




(25) 
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where gi G K[x],hi £ K[x,y]. 

Proof. In the case f\ = h\f[,$2 = ^2/2 having factors in K[x] but gcd(/i,/2) = 1, where 
hi = Cont(/j, y), i = 1, 2, we have 

M(/i, / 2 ) = M(hx, ft) + M(/(, h 2 ) + M(/(, (26) 

Since gcd(/ il /{, /i 2 /*) = 1, M(/ii, /i 2 ) = 0. 

By (|26p . we can assume that Cont(/j,y) = l,i = 1,2. From Lemma 14.11 we have (|2ip . 
We put the triangular sets on the righthand side of (f2T|) into two sets: W\ = M(/fc+i, fk+2) + 
T, k i=1 M(p h f i+1 ), W 2 = EtiMK,/ i+ i) +M(m fc ,/fc+i). It is clear that W 2 C W x . Our 
aim is to delete the multiplicative-zeros of W2 from W\ to derive a group of triangular sets 
with positive coefficients. 

Take any triangular set U = (ui(x),U2(x,y)) out of W2, we can compute gcd of U\(x) 
and some V\(x), where V = {v\{x),v 2 (x,y)) is one triangular set of W\. Denote the gcd as 
p(x). Then decomposing U, we have (p(x),u 2 (x,y)), and put ( ^ff , U2(x, y)) into W2 again. 

Decomposing V, we have (p(x),v 2 (x, y)). And put Cp(x) i V2 i x i n)) ^i- We can compute 
the gcd, say w(x,y), of U2(x,y) and V2(x,y) modulo p(x), then remove the factor w(x,y) 
from u 2 {x,y) (v2(x,y)) and put the left part into ?7(V). We can also use the method in [21] 
to decompose (p(x), U2(x, y)) and (p(x), w 2 (x, y)) into irreducible and regular triangular sets, 
and then we can easily decide that whether two irreducible and regular triangular sets have 
same zero set or not. Thus we can remove the triangular sets in W2. In the end, we can 
remove all the zero sets in W2- We prove the theorem. | 

We will give a multiplicity preserving algorithm to decompose a bivariate polynomial 
system into triangular sets based on the theory above. 

Algorithm 4.7 Input: a zero- dimensional bivariate polynomial system 8?\ = {fi(x,y),f2(x,y)} G 
K[x,y], and d\ = deg(/i,y) > c?2 = deg(/2,y). Output: a group of triangular sets & = 
{[ gi {x),hi{x,y)],i = l,..,n} such that M{f 1 J 2 ) = ^f =1 M{f i ,h i ). 

1. U p = 0,M n = 0. 

2. Compute h { = ContCft, y), % = 1,2. Let /< = fi/h t . M p = U p U {[hi, f 2 ], [Wi]}- 

3. Let mo = 1. While deg(/2,y) > 0, do 

• By pseudo-division, we have mi /1 + gi/ 2 = /3 and /i = Cont(/3, x). f3 = = 
gcd (mi., /i), mi = ^-,ft = J. Let g = gcd(m , h), w = p = ^. If p is not a 
constant, M p = M p U {[p, / 2 ]}. If u> is not a constant, M n = M n U {[w, / 2 ]}. 

• /l = /2, /2 = /3, mo = mi. 
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If deg(/i, y) > 1, M„ = M„ U {[mo, fi}}. 

4. If M n is not empty, following the method in the proof of Theorem 14,61 we can remove 
the multiplicative varieties in IVE n from MLp. Thus W6 obtain a group of triangular sets 
as flZgP - 

Proof. The termination of the algorithm is clear since the degree of fy and fy is finite. The 
correctness of the algorithm is guaranteed by Theorem 14.61 and Lemma 14.11 

Example 4.8 Let C be the curve defined by 

f = 2y 4 - 3y 2 x + x 2 - 2x 3 + x A . 

We will compute the y-critical points (f = ^ = 0) of C. 

Jy = -s- = »y° - 6yx. 
dy 

j- 

Delete the content 2, f y = In the following, we will solve the system E = {f,f y }- 
Following our algorithm, we have 

mif + q 1 f y = pi fy, 
where m\ = 4, p\ = x, q\ = —y, fy = — 3 y 2 + 2 x — 4 x 2 + 2 x 3 . 

m 2 f y + q2fy = P2 h, 

where m,2 = 3, q 2 = ^y,P2 = % (—1 — 16x + 8x 2 ), fy = y. Note that here w 2 = 4. We ignore 
it since it is a constant. 

m 3 fy + q3fy = fy, 

where = 1, q 3 = Sy, fy = 2x{x — l) 2 . Similarly, we ignore u>3 = 3. With Corollary \Jj-2\ 
we have 

M(/, /„) = M(/ 5 , fy) + M(pi,/ W ) + M(p2, fy). 

We can find that 

M(fy,fy) = M(x,y) + m(x-l,y), 

M(pxJ v ) = M(x,4y 3 -3yx) =3M(x,y), 

U(p 2 ,fy) = 2M(x,y)+M(-l-16x + 8x 2 ,-4y 2 + 3x). 

We find that M(x, y) and M(x — 1, y) are zeros with multiplicities 6, 2, respectively. And the 
other zeros are with multiplicities 1. 
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degree 


[5,4] 


[7, 5] 


[9, 7] 


[13. 11] 


[23, 21] 


[33,31] 


MPTD 


0.006 


0.019 


0.105 


1.363 


62.894 


884.577 


BMT 


0.014 


0.024 


0.036 


0.077 


0.280 


0.848 


RC 


0.194 


0.343 


0.810 


4.520 


127.194 


1075.653 


CS 


0.082 


0.880 


25.101 








WS 


0.134 


3.881 


410.399 









Table 1: Timings for different methods 
4.2 Implementation and Comparison 

We implement our algorithm in Maple. We compare the computing time of our methods 
with several other related methods. One is the regular chains method [15] (package 
RegularChains in Maple 13, including two functions), one is Characteristic set method in 
Epsilon[23]. the other is a package wsolve (see [26]). All the results are collected on a PC 
with a 3.2GHz CPU, 2.00G memory, and running Microsoft Windows XP. We use Maple 13 
in the experiments. 

We run 100 examples in each case and compute their average computing time in Table 1. 
We take random dense polynomials with coefficients bounded by [—100, 100] for each exam- 
ple. Table 1 is the timings for the given methods in seconds. Here MPTD means the method 
provided in this paper, RC (BMT) means regular chains method of function "Triangular- 
ize" ( "BivariateModularTriangularize" ) , CS means characteristic set method ("charsets") in 
Epsilon, and WS means wsolve method. The first row is [deg(/i), deg^)]- The first column 
represents the methods. "-" means out of memory or we did not test it. 

We need to mention that only MPTD can compute the multiplicities of the zeros of the 
bivariate polynomial system. BMT and RC are implemented by C code but other's are in 
Maple. MPTD and BMT are only for bivariate polynomial system but other methods work 
for general system. Note that with a mirror modification, MPTD can multiplicity preserving 
triangular decompose system with two multivariate polynomials. 

We can conclude from the table that MPTD and BMT are always faster than RC, CS 
and WS. For the system with low degrees, CS, WS are a little faster than RC, but they both 
are very slow when the degree of the system more than 10. MPTD is always a little faster 
than RC. MPTD is a little faster than BMT for low degree system, but much slower than 
BMT for systems with high degrees. There are several reasons: BMT is in C, using modular 
method and using FFT based arithmetic, but MPTD does not use these techniques. 

5 Conclusion 

We present an algorithm to decompose a polynomial system with two polynomials into tri- 
angular sets. Different from the existing methods for triangular decomposition, our method 
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preserves the multiplicity of the zeros or components of the systems. We implement the 
method for bivariate polynomial systems. We will extend the method to the systems with 
more polynomials in the future. 
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